INVERSE AC JOSEPHSON EFFECT FOR A FLUXON 
IN A LONG MODULATED JUNCTION 



Giovanni Filatrella 

Physikalisches Institut, Universitdt Tubingen 
D-72076 Tubingen, Germany 
filatrella@pit.physik.uni-tuebingen.de 



Boris A. Malomed 
Department of Applied Mathematics, School of Mathematical Sciences 
Raymond and Beverly Sackler Faculty of Exact Sciences 
Tel Aviv University 
Ramat Aviv 69978, Israel 
malomed@math.tau.ac.il 



Robert D. Parmentier 
Dipartimento di Fisica, Universita di Salerno 
I- 8408 1 Baromssi (SA), Italy 
parment @ salerno . inf n . it 



1 



ABSTRACT 



We analyze motion of a fluxon in a weakly damped ac-driven long Josephson junc- 
tion with a periodically modulated maximum Josephson current density. We demonstrate 
both analytically and numerically that a pure ac bias current can drive the fluxon at a 
resonant mean velocity determined by the driving frequency and the spatial period of 
the modulation, provided that the drive amplitude exceeds a certain threshold value. In 
the range of strongly "relativistic" mean velocities, the agreement between results of a 
numerical solution of the effective (ODE) fluxon equation of motion and analytical results 
obtained by means of the harmonic-balance analysis is fairly good; morever, a prelimi- 
nary PDE result tends to confirm the validity of the collective-coordinate (PDE-ODE) 
reduction. At nonrelativistic mean velocities, the basin of attraction, in position-velocity 
space, for phase-locked solutions becomes progressively smaller as the mean velocity is 
decreased. 
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1 INTRODUCTION 



The 'inverse ac Josephson effect' is the generation of dc voltage across the Josephson 
junction by ac bias current flowing through it. This effect was studied in detail for a 
point-like junction 0], as we ^ as f° r shuttle motion of a single fluxon (Josephson vortex) 
in a long but finite junction with reflecting edges |1. In the present work, we consider 
a different dynamical regime which may also be interpreted as the inverse ac Josephson 
effect: steady progressive motion of a fluxon in the ac-driven weakly damped long junction 
with periodically modulated local parameters. In an experiment, this can be realized by 
periodically changing the junction's barrier thickness ||; in the simplest theoretical model 
it is usually assumed that the only local parameter which is modulated is the maximum 
Josephson current, while all others are kept constant ||. Thus, the ac-driven periodically 
modulated junction is described by the normalized, modified sine-Gordon equation in the 
standard notation: 

4>tt — (fixx + sin (ft = e sin(fcx) sin (ft — atfit + 7 sin(atf) . (1) 

where a is the dissipative constant, e and k are the amplitude and wave number of the 
modulation, and 7 and to are the amplitude and frequency of the ac bias current driving the 
junction. We assume that the ac bias current is distributed uniformly along the junction, 
which is a reasonably realistic assumption for overlap- or annular-geometry junctions ||. 

In this work, our objective is to analyze the motion of a fluxon in the model of Eq. (1). 
In the absence of the perturbing factors (e = a = 7 = 0), the fluxon is described by the 
well-known solution of the sine-Gordon equation: 

4>(x, t) = 4 tan 

where V is the fluxon's velocity, and £ = Vt is the coordinate of its center. Following the 
perturbative approach which has proved to be very useful in the theory of long Josephson 
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junctions |7]], we will treat all the terms on the right-hand side of Eq. (1) as small pertur- 
bations. Then, using energy balance or, equivalently (in this case), momentum balance, it 
is straightforward to derive, at the lowest order of the perturbation theory, the following 
equation of motion for the fluxon: 

| f -J==) = ~ V^cos^) - -^L= + sin(^) , (3) 

where, in order to simplify the equation, it is assumed that 



(4) 

i.e., a characteristic width of the fluxon is small in comparison with the modulation 
wavelength (in the opposite limit the coefficient in front of the first term on the right- 
hand side of Eq. (3) would be exponentially small). 

Eq. (3) is an effective equation of motion for a relativistic particle in the presence of 
potential, friction, and driving forces, corresponding, respectively, to the three terms on 
the right-hand side. Since the driving force has no dc (constant) part, it may seem that 
the particle will only oscillate with zero mean velocity. However, in this work we will 
demonstrate that, in the presence of the spatially periodic potential, the ac drive may 
support mean progressive motion of the particle, compensating the dissipation-induced 
braking. Earlier, it has been predicted analytically for models of the Toda-lattice (TL) 
H and of the Frenkel-Kontorova (FK) || type, and shown numerically for the TL model 



in the so-called dual formulation [10| that ac drive may support progressive motion of 
a soliton in weakly damped nonlinear dynamical lattices; a crucial circumstance which 
made the effect possible was the spatial periodicity of the lattice. Moreover, the stable 
propagation of ac-driven solitons in a model electric transmission line described by the 
damped dual TL equation was detected experimentally in the work [ 11]] . Note that for the 
TL in its usual form the effect does not take place because the ac drive can provide only 



for the balance of energy, but not for that of the momentum, while in the dual formulation 



the soliton's momentum is identically zero, and energy balance becomes sufficient |TI|. In 
the FK model the balance of energy alone is, in principle, sufficient since the momentum 
is not conserved in this system even in the absence of dissipation and drive, but the 
threshold predicted by the theory M for the effect proves to be so high that it would be 



hard to observe in a numerical experiment [12 



In the present work, we are dealing with the continuum model of Eq. (1). However, 
in the presence of a periodic potential, the ac drive allows the soliton in this damped 
continuum system to move at a nonzero mean velocity, the underlying mechanism being 
essentially the same as in the discrete lattices. Indeed, if the soliton is moving with a 
mean velocity V through the periodic potential, this gives rise to the temporal frequency 
kV (recall that k is the wave number which determines the spatial periodicity). The 
condition 

kV = u , (5) 

where u is the driving frequency, may give rise to a resonant transfer of energy from the 
ac driving field to the soliton. In turn, this supply of energy gives the soliton a chance to 
compensate the dissipative losses. In the case of the long Josephson junction, the mean 
dc voltage across the junction is proportional to the mean velocity of the fluxon. Thus, 
supporting the fluxon's progressive motion by the ac bias current may be regarded as 
another manifestation of the inverse ac Josephson effect. 

In section 2 of this work, we will develop a fully analytical approach to investigate the 
ac-driven motion of the fluxon in the "ultrarelativistic" case, i.e., when the mean velocity 
is close enough to the junction's limit (Swihart) velocity (which is equal to one in the 
notation adopted). In section 3, we will display results of the full numerical solution of 
the ODE, Eq. (3). 



2 ANALYTICAL CONSIDERATIONS 

An analytical approach to Eq. (3) is possible if one assumes that the small parameters 
a and 7 are "very small", while e is "not so small" (the physically relevant range for e 
is < e < 1). In other words, we will consider, in the lowest approximation, motion 
of the free relativistic particle in the periodic potential, and then the friction and driving 
forces will be taken into account as additional small perturbations. Another fundamental 
assumption which makes the analytical consideration possible is that the particle is moving 
at ultrarelativistic velocities, i.e., we will set 

i = l-w 2 ,w 2 <1. (6) 

Insertion of Eq. (6) into Eq. (3) leads, after simple algebra, to the following reduced 
equation of motion for the variable w: 

^ = ekw 3 cos(£;£) + aw ^= w 2 sin(ujt) . (7) 

at 2y 2 

In what follows below, we will concentrate on the case when the "particle" moves with 
a nearly constant velocity. In this case, the variable £ in the argument of the cosine in the 
first term on the right-hand side of Eq. (7) may be replaced, in the lowest approximation, 
by 

e (0) = V t + £ , (8) 

where V is the mean velocity (the same as in Eq. (5)), and £ is a constant phase shift 
between the particle's mean law of motion and the time-periodic driving force. Finally, 
introducing the renormalized time r = kVt, we obtain the lowest-order equation of motion 
in the form (with 8 = fc£o) 

dw e n a 777 2 . u . 

— = = w 6 cos(t + 8) + w sin . (9) 

dr V K J kV 2^2k K kV ' y ' 



Now, we will make use of the fundamental assumption mentioned above, viz., that 
the driving and friction terms may be regarded as small in comparison with the first 
(potential) term on the right-hand side of Eq. (9). In the zero'th-order approximation, 
we drop the small terms, obtaining from Eq. (9) 

^ = ^w 3 cos(r + 5). (10) 
ar V 

The exact solution of Eq. (10) is 



w = 



A 

\j 1 - M sin (r + 5) ' (U) 

where A is an arbitrary constant. According to Eq. (6), w 2 <C 1, which implies that 
A <C 1. Moreover, the physically significant range for e is e < 1; and finally, V ~ 1. 
Consequently, the coefficient of the sine term in the denominator of Eq. (11) is small 
enough to replace this expression by its expansion, 



w 



1 + ■=■ sm(r + 5) 



(12) 



Eqs. (6) and (12) describe the motion of a particle at the mean velocity V = 1 — A, 
with superimposed small oscillations of the velocity corresponding to the second term in 
the expression (12). Inserting this V into the resonance condition of Eq. (5), we find the 
following relation between the driving frequency uj and the corresponding resonant value 
of the parameter A: 

u = k(l-A). (13) 

Obviously, in a real physical situation the system chooses itself the resonant value of A 
corresponding to given u. However, in the analytical consideration that follows below it 
will be more convenient to regard A as a given arbitrary parameter (small enough), and 
then to match the frequency to it according to Eq. (13). 
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A resonant drive can support progressive motion of a particle in a lossy periodic 
medium if the drive amplitude exceeds a certain minimum (threshold) value proportional 
to the friction coefficient ||. In fact, the threshold value of the amplitude, regarded as 
a function of the driving frequency, is the basic characteristic of such ac-driven motion 



of the soliton in the damped system || 1 1 . In the present case, the amplitude of the 



driving force is the parameter 7. A convenient, approximate technique for calculating 



its threshold value, 7thr, is to use the idea of harmonic balance |13| , i.e., we insert the 
approximate solution given by Eq. (12) into the last two terms on the right-hand-side of 
Eq. (9), adjusting the coefficients in these terms so that the zero'th harmonic (and, in 
general, as many harmonics as possible) vanishes. The result of such a calculation (for 
the zero'th harmonic) is 

7 = l^ A -3/2 / | cos(J) |_ (M) 

Recalling that 5 is an arbitarary phase shift, it is apparent that the threshold corresponds 
to I cos(5)| = 1 in Eq. (14), i.e., 

7thr = — A" 3 / 2 . (15) 

Tie 

This is our main analytical result. The proportionality of 7 t h r to the dissipation 
constant a, as well as the inverse proportionality to e, are obvious. A nonobvious feature 
of Eq. (15) is the power — | of A (it is evident, of course, that this power must be 
negative). In the next section, comparing with results of direct numerical simulations of 
the equation of motion, Eq. (3), will show that this particular value of the power, as well 
as the numerical prefactor in Eq. (15), are quite reasonably accurate. 
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3 NUMERICAL RESULTS 

In the numerical integration of Eq. (3), we fixed the driving frequency to ou — 1. The 
justification for this step is provided by noting that u> can be scaled out of Eq. (3), 
i.e., set to unity, by defining the following re-scaled variables: t' = ut, £' = (which 
leaves £ unchanged), k' = k/u, a' = a/uj, 7' = 7/0;, and e ; = e. Such a re-scaling is 
valid provided that condition (4) continues to hold, but it obviously neglects any possible 
frequency-dependent resonance effects with the background oscillation (for the case when 
the driving frequency is close to the junction's plasma frequency, the resonance effects were 



analyzed in detail in |14|)- Then, with A taken as an arbitrary parameter (small enough), 
we calculated the modulation wave number k from the resonant relation (13). The mode 
of simulations chosen was as follows: with all the parameters but the drive amplitude 7 
fixed, we reduced the value of 7 quasi- adiabatically until the solution dropped out of 
synchronism with the driver; the smallest value of 7 so obtained was taken as 7thr • The 
numerical results obtained are summarized in Figs. 1 and 2. 

Figs. 1(a) and 1(b) show, respectively, the inverse and direct proportionality of 7 t h r 
to the parameters e and a. The analytical prediction of Eq. (15) is qualitatively well 
confirmed, even though there is a small quantitative discrepancy between analytical and 
numerical results. Fig. 2, instead, shows the numerical dependence on A. Here, we see 
that, for small enough values of A, the agreement between analytical and numerical results 
is quite good. The agreement gradually worsens as A increases, which is not surprising in 
view of the fact that A <C 1 is a fundamental assumption of the ultrarelativistic analysis 
leading to Eq. (15). 

In the extreme nonrelativistic limit, i.e., A — > 1, which implies V — > 0, Eq. (3) 
becomes, with trivial modifications, the equation for the damped, ac-driven, simple pen- 
dulum. The existence of a threshold value of the drive amplitude for phase-locked motion 
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in this case was recently derived by Filatrella et al. @] using an energy-balance argument. 
The solid curve in Fig. 2 shows the drive-amplitude threshold calculated for this case 5 as 
is evident from the figure, the numerical solution of Eq. (3) increasingly approaches this 
curve as A — > 1. 

In Ref. |J, the numerical solutions of the "nonrelativistic" equation showed a gradual 
shrinkage, and eventual disappearance, with decreasing u, of the basin of attraction in 
position-velocity space of the characteristic phase-locked motion. We have seen the same 
phenomenon, which, from Eq. (13), corresponds to the situation A — > 1, in our present 
numerical solutions of Eq. (3). We suggest that this phenomenon might be attributable to 
the co-existence, in this region of parameter space, of chaotic trajectories in the underlying 
dynamics. 

A crucial step in our analysis has been the reduction of the system description from 
the PDE of Eq. (1) to the ODE of Eq. (3). In order to check that the existence of 
ac-driven phase-locked states is not just an artifact of this reduction we have performed 
a very preliminary numerical study of Eq. (1), in which we have seen clear evidence, 
for several different parameter combinations, that the phenomenon is also present in the 



PDE In these preliminary tests we applied periodic boundary conditions to Eq. (1), 
with the junction length chosen to contain one period of the spatial modulation, i.e., 
L = 2ir/k, with k given by Eq. (13). The PDE was discretized in space using a three- 
point approximation for the second spatial derivative, and integrated in time using a 
fourth-order Runge-Kutta routine. Eq. (2), with £ = and V = 0, was used as the initial 
condition. In one such run, with a = 0.005, e = 1.0, 7 = 0.53, uj = 1.0, A = 0.1 and with 
a time step of 0.02 and a spatial grid of 100 points, the fluxon continued to propagate for 
at least 2500 normalized time units. Further detailed support for the very broad range of 
validity of the collective-coordinate (PDE-ODE) reduction in the periodically modulated 
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sine-Gordon equation (without loss or driving terms) has recently been presented by 
Sanchez et al. [|16 . 



4 CONCLUSIONS 

In this work, we have demonstrated that a soliton in a periodically inhomogeneous con- 
tinuum lossy medium can be driven at a nonzero mean velocity by an ac driving force. 
We have also found a good agreement between the analytical approximation based on 
harmonic balance and direct numerical simulations of the effective (ODE) equation of 
motion for the soliton. A noteworthy result is that the soliton can be easily driven in 
the "ultrarelativistic" region, but less easily at moderately relativistic and nonrelativistic 
values of the mean velocity. In this connection, one might reasonably ask, if the basin 
of attraction for ac-driven phase-locked states is small, whether or not such states can 
be of any potential physical interest, e.g., for practical applications. Our answer to this 
question is positive: the fact that an experimentalist might have to "prepare" his system, 
e.g., by initially applying a combination dc + ac drive, and then gradually reducing the 
dc bias to zero, is not a serious impediment. What is important is that the basin of 
attraction be sufficiently large so that, once having arrived in an ac-driven phase- locked 
state (no matter how), the system be sufficiently stable against thermal noise and other 
perturbations, for a sufficiently long time to make measurements. The technique of cell- 
to-cell mapping |[7|, [181 , applied to Eq. (3), might provide a useful tool for obtaining some 



initial insight into this problem. 

Finally, it is relevant to mention that in this work we have confined ourselves to 
consideration of the simplest resonance only. Higher-order resonances may also take place, 
in which mkV = nuo, with incommensurable integers m and n, cf. Eq. (5). It is natural 
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to expect that the threshold amplitudes for the higher resonances will be considerably 
larger than for the fundamental resonance considered here. 
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Figure Captions 

Fig. 1. Dependence of the threshold value of the driving force, 7 t h r , on the parameters: 
(a) e (modulation coefficient), with a = 0.001, oo = 1.0, A = 0.1; (b) a (dissipation 
coefficient), with e = 1.0, u = 1.0, A = 0.1. Dash-dotted lines: prediction of Eq. (15); 
asterisks: numerical solutions of Eq. (3). Solid lines connecting the asterisks are a guide 
to the eye only. 

Fig. 2. Dependence of the threshold value of the driving force, 7 t hr, on the velocity pa- 
rameter, A, with a = 0.001, uj = 1.0, e = 0.5. Dash-dotted line: prediction of Eq. (15); 
asterisks: numerical solution of Eq. (3); solid curve: analytical prediction of a nonrel- 
ativistic theory (see text). Solid lines connecting the asterisks are a guide to the eye 
only. 
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